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Abstract. We present a multivariate generating function for all n X n non- 
negative integral matrices with all row and column sums equal to a positive 
integer t, the so called semi-magic squares. As a consequence we obtain for- 
mulas for all coefficients of the Ehrhart polynomial of the polytope B n of n X n 
doubly-stochastic matrices, also known as the Birkhoff polytope. In particular 
we derive formulas for the volumes of B n and any of its faces. 



1. Introduction 

Let B n denote the convex polytope of n x n doubly-stochastic matrices; that 
is, the set of real nonnegative matrices with all row and column sums equal to 
one. The polytope B n is often called the Birkhoff-von Neumann polytope, the 
assignment polytope, or simply the Birkhoff polytope. It is a well-known problem 
to compute the volume of B n and there is a fair amount of work on the topic 
(see [5j [TTJ [15] and the references therein for information on prior work); in this 
paper, we present the first exact formula for the volume of B n . The formula will 
follow from a multivariate rational generating function for all possible n x n integer 
nonnegative matrices with all row and column sums equal to a positive integer t, 
the so called semi-magic squares [161 122] (although many authors refer to them as 
magic squares). 

Before stating our main formula, we give a few necessary definitions and notation. 
We call a directed spanning tree with all edges pointing away from a root i an l- 
arborescence. The set of all £-arborescences on the nodes [n] = {1, 2, . . . , n} will be 
denoted by Arb(^,n). It is well known that the cardinality of Arb(£, n) is n n ~ 2 . 
For any T £ Arb(£, n), we denote by E(T) the set of directed edges of T. As usual 
let S n be the set of all permutations on [n]. For any a € S n , wc associate a with 
its corresponding permutation matrix, i.e., the n x n matrix whose (i,o~(i)) entry 
is 1 and zero otherwise. Throughout this paper, we will use a to denote both a 
permutation and the corresponding matrix and it should be clear which one it refers 
to according to the context. The bracket operator (•, •) denotes the dot product of 
two vectors. 

It is well known that given a <i-dimensional integral polytope P, that is a polytope 
whose vertices have integer coordinates, for any positive integer t, the number e(P, t) 
of lattice points contained in the i-th dilation, tP = {tX | X e P}, is a polynomial 
of degree dim(P) in the variable t. Furthermore, the leading coefficient of e(P, f) 
is the normalized volume of P in units equal to the volume of the fundamental 
domain of the affine lattice spanned by P (see Chapter 4 of [22] or the book [7]). 
This polynomial is called the Ehrhart polynomial of P. 

l 
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One can find an expression for the Ehrhart polynomial e(B n ,t) of B n using the 
multivariate generating function 

f(tB n ,z) = zM 

MetB n nz™ 2 

of the lattice points of tB n , where z M = Yii<i j< n z Tj ° ^ M = (TOij) is an n by n 
matrix in M n . One can see that by plugging = 1 for all i and j in f(tB n ,z), 
we get e(B n ,t). Our main result is 

Theorem 1.1. Given any positive integer t, the multivariate generating function 
for the lattice points of tB n is given by the expression 

a- 1 ) /(«*».»)=£ E II n-uz^y 

f£S„ TeArb(£,n) e£E(T) y 11 



"k,a(k) ■ 

Here W T ' e denotes the n x n (0, —1, l)-matrix associated to the unique oriented 
cycle in the graph T + e (see Definition \ 3.17\ for details) and W T ' e o~ denotes the 
usual matrix multiplication ofW T ' e and the permutation matrix a. 

As we apply Lemma T5. 41 to Theorem ll.il we obtain the desired corollary: 

Corollary 1.2. For any choice of fixed I £ [n], the coefficient of t k in the Ehrhart 
polynomial e(B n ,t) of the polytope B n of n x n doubly- stochastic matrices is given 
by the foi~mula 

n ox I V V (^^)) k td d - k ({{c,W T ' e a), e£E{T)}) 

ceS n TeArb(l,n) llefh(I)\ ' / 

In the formula W T,e is the n x n (0, —1, l)-matrix associated to the unique oriented 
cycle in T + e as defined in Definition \3.17\ and W T,e o~ denotes the usual matrix 
multiplication of W T ' e and the permutation matrix a. The symbol tdj(S') is the 
j-th Todd polynomial evaluated at the numbers in the set S (see Definition \5.1\ for 
details). Finally, c £ K™ is any vector such that (c,W T ' e a) is non-zero for all 
pairs (T, e) of an I- arborescence T and a directed edge e ^ E(T) and all a £ S n . 
As a special case, the normalized volume of B n is given by 

i / \(n-l) 2 

(L3) vol(B„) = jj—Jyry E E n (cWT^y 

We stress that each rational function summand of Formula (|1.1|) is given only in 
terms of trees and cycles of a directed complete graph. Our proof of Theorem 11.11 
is based in the lattice point rational functions as developed in [4] with some help 
from the theory of Grobner bases of toric ideals as outlined in [23) . 

There is a large collection of prior work on this topic that we mention now to 
put our result in perspective. In [S] the authors computed the exact value of the 
volume and the Ehrhart polynomials for up to n = 10, which is the current record 
for exact computation. The computations in [5] took several years of computer CPU 
(running in a parallel machine setup) and our volume formula is so far unable to 
beat their record without a much more sophisticated implementation. On the other 
hand, in two recent papers, Canficld and McKay [9j [TO] provide simple asymptotic 
formulas for the volume of B n as well as the number of lattice points of tB n . 
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However, our closed formula for the volume of B n is nonetheless interesting for 
the following reasons. First, as it was demonstrated in [HQ], the faces of B n are 
also quite interesting for combinatorics and applications. For example all network 
polytopes appear as faces of a large enough B n . From our formula it is easy to work 
out volume formulas for any concrete face of B n . We demonstrate this possibility 
in the case of the well-known CRY n polytope [12] whose volume is equal to the 
product of the first n — 1 Catalan numbers (see [25] ). Concretely we obtain for 
the first time the Ehrhart polynomials of facets of B n and CRY n for n < 7. In 
principle, this could be applied to derive formulas for the number of integral flows 
on networks. Second, not only we can derive formulas for the coefficients of the 
Ehrhart polynomial of B n , but we can also derive formulas for the integral of any 
polynomial function over B n . We hope our generating function will be useful for 
various problems over the set of all semi-magic squares, at least for small values of 
n. 

This paper is organized as follows: In Section [2] we begin with background ma- 
terial that will be used in forthcoming sections, including background properties of 
B n , a short discussion of Grobner bases and triangulations, Brion's theorem and 
generating functions for lattice points in polyhcdra. In that section, we sketch 
the steps we will follow to compute the generating function of lattice points inside 
cones. In Section [3] we discuss the triangulations of the dual cone at each vertex of 
B n which we encode via Grobner bases. From Brion's formula we derive in Section 
2] a sum of rational functions encoding all the lattice points of the dilation tB n 
and thus a proof of Theorem 11.11 In Section [5] we show how from Theorem 11.11 we 
can derive all the coefficients of the Ehrhart polynomial of B n after expressing the 
generating function in terms of Todd polynomials. Finally, in Section^ we explain 
how to obtain Ehrhart polynomials and formulas of integration for any face of B n . 

2. Background 

For basic definitions about convex polytopes which are not stated in this paper, 
please see [26,. Chapters 5 and 6 in [24] have a very detailed introduction to B n 
and transportation polytopes. For all the details and proofs about lattice point 
counting and their multivariate generating functions see [31 HI E] . We begin with 
some useful facts about the polytope B n . It is well known that the vertices of B n are 
precisely the n x n permutation matrices. Permutation matrices are in bijection 
with matchings on the complete bipartite graph K n ^ n . The polytope B n lies in 
the n 2 -dimensional real space I™ = {n x n real matrices}, and we use M(i,j) 
to denote the (i,j)-entry of a matrix M in the space. There is a graph theoretic 
description of the edges of B n ; they correspond to the cycles in K n ^ n . On the other 
hand, for each pair with 1 < i, j < n, the set of doubly-stochastic matrices 
with (i,j) entry equal to is a facet (a maximal proper face) of B n and all facets 
arise in this way. It is also easy to see that the dimension of B n is (n — l) 2 (i.e., 
the volume we wish to compute is the (n — l) 2 -volume of B n regarded as a subset 
of n 2 -dimcnsional Euclidean space). Note that an n x n doubly-stochastic matrix 
is uniquely determined by its upper left (n — 1) x (n — 1) submatrix. The set 
of (n — 1) x (n — 1) matrices obtained this way is the set A n of all nonnegativc 
(n — 1) x (n — 1) matrices with row and column sums < 1 such that the sum of all 
the entries is at least n — 2. A n is afhnely isomorphic to B n and we often compute 
in A n instead of B n because A n is full-dimensional. 
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Cones and Generating functions for lattice points. For any polytope P £ R , 
we would like to write a generating function for the following sum encoding the 
lattice points of P 



where z Q = z" 1 z^ 2 ■ ■ ■ z% d . We give now a step- by-step description of how the 
generating function is constructed. 

A cone is the set of all linear nonncgative combinations of a finite set of vectors. If 
a cone contains no other linear subspace besides the origin then we say it is pointed. 
Given a cone C C R d , the dual cone to C is a cone C* = {y E R d | (x, y) > 0, Mx £ 
C}. The following lemma states some properties of dual cones. (See Theorem 9.1 
in [21] for a proof). 

Lemma 2.1. Let C be a pointed cone in R™, and let D = C* be its dual cone. 
Then the following properties hold 

(1) C is the dual cone of D, namely C = D* = (C*)*. 

(2) If C is a full dimensional pointed cone, then so is D. Moreover, if {Fi} is 
the set of facets of C, then D is precisely the cone generated by the set of 
rays {Ri} satisfying, for any i, 

(2.1) Ri is perpendicular to Fi, and for any ray R of C not on Fi : (R,Ri) > 0. 

Now if P is a polytope and v is a vertex of P, the supporting polyhedron of P at 
v is 

S(P, v) = v + {u £ R. d : v + Su £ P for all sufficiently small S > 0}, 

and the supporting cone of P at v is defined as C{P, v) = S(P, v) — v. 
For a set icR 11 , the indicator function [A] : R d — ■> R of A is defined as 



The algebra of polyhedra P(R d ) is the vector space over Q spanned by the indicator 
functions [P] of all polyhedra P C R d . The algebra of polytopcs P P (R d ) is the 
subspace spanned by the indicator functions of the polytopes in R d . The algebra of 
cones Pc(R d ) is the subspace spanned by the indicator functions of the polyhedral 
cones in R rf . A linear transformation 



where V is a vector space over Q is called a valuation. Similarly, linear transfor- 
mations defined on Pp(R d ) and P c (R d ) arc also called valuations [4]. 

One important tool for counting lattice points is the ability of expressing the 
indicator function of a simplicial cone as an integer linear combination of the indi- 
cator functions of unimodular simplicial cones. Given a cone K C R d , we say that 
the finite family of cones Ki, i £ I = {1, 2, . . . , 1} is a decomposition of K if there 
are numbers e; G { — 1,1} such that 




z 



.0 



a£PDZ d 




$ : P(R d ) -> V, 



[K]=J2ti[Ki]. 



Theorem 2.2 (Theorem 3.1 and its proof in [4]). There is a map J which, to 
each rational polyhedron P C R d , associates a unique rational function f(P, z) in 
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d complex variables z 6 C d , z = (zi, . . . , suc/i i/iai i/ie following properties are 
satisfied: 

(i) T/ie map $ is a valuation. 

(ii) J/ P is pointed, there exists a nonempty open subset U p C C d , such that 
J2aePnz d z " converges absolutely to f(P,z) for all z G Up. 

(iii) If P is pointed, then f(P, z) satisfies 

f(p,z)= J2 z " 

aGPnZ d 

/or any z £ C d where the series converges absolutely. 

(iv) If P is not pointed, i.e., P contains a line, then f(P, z) = 0. 

Because the rational function f(P, z) encodes the lattice points of P, we call 
/(P, z) the multivariate generating function of the lattice points (MGF) of P. The 
rational function has an expression as a sum of simple terms, but to describe them 
we need the following facts. 

Theorem 2.3 (Brion, 1988; Lawrence, 1991). (see [H [3] for proofs) Let P be a 
rational polyhedron and let V(P) be the vertex set of P. Then, 

f(P,z)= ]T f(S(P,v),z). 

This theorem reduces the problem of finding the MGF of a rational polyhedron 
P to that of finding the MGF of the supporting polyhcdra at each vertex of P. 
If the vertex of the supporting polyhedron is integral we can simply assume the 
vertex is the origin and work instead with supporting cones. 

Corollary 2.4. If P an integral polyhedron, i.e., all the vertices of P are 
integral vertices, then 

f(P,z)= £ 7?f{C{P,V),z). 

vev(P) 

Although it is in general more complicated to give the MGF of an arbitrary cone, 
if the cone is unimodular, its MGF has a simple form: 

Lemma 2.5 (Lemma 4.1 in [5]). If K is a d-dimensional pointed cone in W l 
generated by the rays {ri}i<i<d, where the Ti 's form a %-basis of the lattice Z" n 
span(K) (span(K) is the d-subspace where K lies), then we say K is a unimodular 
cone and we have that 

Barvinok gave an algorithm to decompose any pointed cone C as a signed sum 
of simple unimodular cones [3] and thus deriving an expression for f(P, z) as a sum 
of terms like those in Lemma 12.51 In principle, one needs to keep track of lower 
dimensional cones in the decomposition for writing a inclusion-exclusion formula of 
the MGF f(C, z). Fortunately, by using the Prion's polarization trick (see Remark 
4.3 in [4]), one only needs to consider full-dimensional cones. This trick involves 
using dual cones of a decomposition of the dual cone to C instead of directly 
decomposing C. The main idea is to note that the duals of low dimensional cones 
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are not pointed and thus, from Part (iv) of Theorem 12. 2\ their associated rational 
functions vanish. 

Now we are ready to sketch the main steps of Barvinok's algorithm to compute 
f(C,z) (sec [J for details): 

Algorithm: 

Input: a rational cone full dimensional pointed cone C. 
Output: the MGF of C : /(C, z). 

(1) Find the dual cone D = C* to C. 

(2) Apply the Barvinok decomposition to D into a set of unimodular cones Di 
which have the same dimension as D (ignoring all the lower dimensional 
cones) . 

(3) Find the dual cone Ci to each Di. The cone Ci will be unimodular as well. 

(4) f(C,z) = Yli e i/(Ci) z )j where is +1 or —1 determined by Barvinok 
decomposition. 

This algorithm is still not right for us; the algorithm is for full-dimensional cones, 
however, the cones we need to study arc not full-dimensional since the Birkhoff 
polytope is not full-dimensional. Also, Lemma 12.11 provides us a way to compute 
the rays of C* if C is full dimensional and pointed. Hence, it will be nice if we can 
make our cones full-dimensional. What we will do is to properly project cones into 
a lower dimensional space so that they become full-dimensional. 

Definition 2.6. Let V C R™ and W C R m be vector spaces with full rank lattices 
Ly := V n Z" and Lw := W (~l Z m , respectively. A linear map </> from V to W is a 
good projection if tp gives a bijection between Ly and Lyy. Note that because of 
the linearity of <p, the lattices Ly and Lw have the same rank. 

Lemma 2.7. Suppose V, W are as in Definition \2.6\ and </> is a good projection from 
V to W. Because <f> is a linear map, we can consider 4> is given by a certain m x n 
matrix (f> — {(f>ij). We define a map $ : C" 1 — > C" by mapping y = (j/i, . . . ,y m ) G 
C m toz, = (z u \ . . , z n ) G C™, where 

m 

Then the following statements hold. 

(1) dim(V) = dim(W). 

(2) 4> gives an isomorphism between V and W which preserves the lattice. 
Therefore, there exists an inverse (linear) map (f>~ 1 from W to V that 
preserves the lattice as well. Thus, 4>~ l is also a good projection from W 
to V. 

(3) C is a unimodular cone in V if and only if 4>{C) is a unimodular cone in 
W. 

(4) For any a e Z n and y G C"\ if = 4>{a) and z = $(y), then y fj = z Q . 

(5) For any pointed rational polyhedron P G V, the series J2pe<p(P)nz m con ~ 
verges absolutely if and only if the series J2 a ePnz n ^(y) a converges abso- 
lutely. Furthermore, we have 



(2.2) 



/(0(P),y)=/(P,$(y)). 
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(6) Let P±, P2, . . . , Pk be pointed rational polyhedra in V, and a\,...,ak £ C, 
then 

f(P,z) =^a i /(P i) z)^/(0(P),y)=^a i /(</»(P i ),y)- 

Proof. The proofs of (1), (2) and (3) follow from the fact that good projections 
give lattices of same rank and thus isomorphic vector spaces. For the proof of (4), 
/3 = 4>{a) implies that Pi = 5Zj=i 4>i,j a j- Thus, 

n n rn m n m m 

- II OH vt« a > 1111 yf^ = 11 «F " jaj = 11 # = y". 

j — l J=l i— 1 i— 1 i=l i=l i— 1 

Because (pis a good projection, the lattice in P and the lattice in 4>(P) are in one- 
to-one correspondence under </>. Therefore, to prove (4), it is enough to show that 
if (3 — 4>{a) and z = &(y), then y 13 = z a . [3 = 4>{a) implies that Pi = 52j=i 4>i,j a 3- 
Thus, 

n n rn m n m rn 

<r II z? 1111 1111 y*-*" = II yf U = H yf = y*. 

j — l j — li—1 i— 1 j — 1 i— 1 i— 1 

The first part of (5) follows immediately from (4). Let Y be the set of y £ C m for 
which the series E/3e0(p)nz m V 13 converges absolutely and Z be the set of z € C™ 
for which the series Eagpnz™ zQ converges absolutely. By the first part of (4), 
C Z. By Theorem /(^z) = E QS p n z" z ° for an y z e Z - In particular, 
/(P z ) = E ae pnz™ z " for an y z e *( y )- Hence, for any y E Y. 

/(p,$( y ))= Yl w = E y' 3 - 

oePnZ" /3e0(P)nz™ 

We use Theorem 12.21 again to conclude that /(P, $(?/)) is the rational function 
/(</>(P),y) associated to </>(P). 

Given (2), we only need to check one direction in (6). Suppose /(P, z) = 
X>i/(Pi,z). We can apply $T2$ on both sides to obtain /(0(P), y) = £ a;/(</>(P), y) 

□ 

Using Lemma |2.7[ wc modify Barvinok's algorithm and sketch a method to con- 
struct /(C, z) for supporting cones C at vertices of B„ . We will try to follow this 
sequence of steps in Section [3J 

(CMGF) Method for constructing the multivariate generating function 
for lattice points of a cone: 

Input: a rational (not necessarily full dimensional) pointed cone C C W 1 . 
Output: the MGF of C : f(C,z). 

(0) Let V be the subspacc spanned by C in 1™. Find a subspace W of W m 
together with a good projection <f> from V to W. Let C = <f>(C). 

(1) Find a dual cone D to C. 

(2) Decompose D into addition and subtraction of unimodular cones Di which 
have the same dimension as D, ignoring all the lower dimensional cones. 

(3) Find dual cone d of each Di. Note, that Ci is also unimodular. Let 

a = <t>-\c l ). 

(4) f(C,z) = ^2 i eif(Ci,z), where is +1 or —1 determined by the signed 
decomposition. 
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In the next section, we will apply the method (CMGF) step by step to the 
supporting cone at the vertex J, the identity permutation. We will get the MGF 
of this supporting cone and, by applying the action of the symmetric group 

£>n ? we 

can deduce the MGF of all other supporting cones of vertices of B n and thus, by 
Theorem 12.31 the MGF of B n . We will see later, in Section [5l that the knowledge 
of /(P, z) as a sum of rational functions yields a rational function formula for the 
volume of P. 

Triangulations and Grobner bases of toric ideals. For step (2) in our step-by- 
step construction of the generating function, we will show (Lemma l3.4[) that in fact 
any triangulation of the dual cone of the supporting cone of a vertex gives already 
a set of unimodular cones (hence, the e^'s in Step (4) are all +1). A triangulation 
of a cone C is a special decomposition of a cone as the union of simplicial cones 
with disjoint interiors whose union covers completely the cone C. In this article 
we use polynomial ideals to codify the triangulations, namely toric ideals and their 
Grobner bases. See Chapter 8 in [23] for all details. Here are the essential notions: 
Fix a set A — {ai, 0,2, . . ■ , a n } C 7L d . For any u = (u-i, 1*2, ■ ■ ■ , u n ) € Z n , we let 

u.4 := mai + u 2 a 2 + h u n a n . 

For any u G Z d , we denote by supp(u) := {i \ m / 0} the support of u. Every 
u e 7L d can be written uniquely as u = u + — u~, where u + and u~ are nonnegative 
and have disjoint support. 

Definition 2.8. The toric ideal of A, Ia C k[x\ := k[x\, X2, ■ ■ ■ ,x n ] is the ideal 
generated by the binomials 

I A := (x u+ -x u ~ I u„4 = 0). 

Given a real vector A = (Ai, . . . , A„) in R™, we can define a monomial order >\ 
that for any a, b e Z> , their monomials satisfy x a >\ x b if (a, A) > (6, A) and ties 
are broken via the lexicographic order. Using the ordering of monomials we can 
select the initial monomial of a polynomial / with respect to >a, i.e., the highest 
term present. We will denote it by in >x (f). For an ideal I contained in C[xi, .., x n ] 
its initial ideal is the ideal in >x (I) generated by the initial monomials of all polyno- 
mials in I. A finite subset of polynomials G — {gi, g n } of an ideal / is a Grobner 
basis of I with respect to >a if in> x (I) is generated by {in >x (gi), in >x (g n )}- In 
other words, G is a Grobner basis for / if the initial monomial of any polynomial in 
/ is divisible by one of the monomials in >x (gi). It can be proved from the definition 
that a Grobner basis is a generating set for the ideal /. As we will state later, each 
Grobner basis of the toric ideal Ia yields a regular triangulation of the convex hull 
of A. The fact that triangulations constructed using Grobner bases are regular will 
not be used in our construction. 

A subdivision of A is a collection T of subsets of A, called cells, whose convex 
hulls form a polyhedral complex with support Q = conv(A). If each cell in T is 
a simplex, then T is called a triangulation of A. Every vector A = (Ai, . . . , A„) in 
W l induces a subdivision of A = {ai, . . . ,a n } as follows. Consider the polytope 
Qx = conv({(ai, Ai), . . . , (a n , A„)}) which lies in R d+1 . Generally, Q\ is a polytope 
of dimension dim(conv(A)) + 1. The lower envelope of Q\ is the collection of faces 
of the form {x S Q\\(c, x) = cq} with Q\ contained in the halfspace {c,x) < co 
and the last coordinate Cd+i is negative. The lower envelope of Q\ is a polyhe- 
dral complex of dimension dim(conv(A)) . We define T\ as the subdivision of A 
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whose cells are the projections of the cells of the lower envelope of Q\. In other 
words, {a ii; Oj 2 , . . . , a ik } is a cell of T\ if {(a^, A,J, (o, 2 , Aj 2 ),..., (a, fe , AiJ} arc 
the vertices of a face in the lower envelope of Qa- The subdivision T\ is called 
a regular subdivision of .A. Remark that just as a triangulation can be uniquely 
specified by its maximal dimensional simplices, it can also be uniquely expressed 
by its minimal non- faces (minimal under containment). Now we are ready to state 
the algebra-triangulation correspondence: 

Theorem 2.9 (See proof in Chapter 8 of [13]). Let A be an n x d matrix with 
integer entries, whose rows vectors {a±, . . . , a n } span an affine space of dimension 
d — 1. Let I _4 be the toric ideal defined by A. Then, the minimal non-faces of the 
regular triangulation of A associated to the vector A can be read from the generators 
of the radical of the initial ideal of the Grobner basis of I a with respect to the term 
order >\. More precisely, for A generic, the radical of the initial ideal of 1 a equals 

: {ii,«2, ■ • - ,i s } is a minimal non-face of T\) = Q (xi : i £ a). 

The crucial fact we will use is that the maximal simplices of the regular trian- 
gulation T\ are transversals to the supports of the monomials from the initial ideal 
of the Grobner basis. In the next section, we will apply Theorem 12.91 to create a 
triangulation of the dual cones. 

To the readers who are unfamiliar with commutative algebra language, using a 
Grobner basis to describe a triangulation may not feel totally necessary or clear. 
Thus, we explain here the advantages of doing it this way. First, traditionally 
checking that a set of simplices is a triangulation of A is not trivial since one has 
to verify they have disjoint interiors (which requires a full description of all linear 
dependences of the rays) and that the union of the simplicial cones fully covers the 
convex hull of A. But, having a Grobner basis avoids checking these two tedious 
geometric facts. Second, the initial monomials of the Grobner bases are precisely 
the minimal non- faces of the triangulation T\, which are complementary to the 
maximal simplicial cones of the triangulation. From the point of view of efficiency, 
the encoding of a simplicial complex via its non-faces is sometimes much more 
economic than via its maximal facets. For more on the theory of triangulations see 

3. The MGF of the supporting cone of B n at the vertex / 

Due to the transitive action of the symmetric group on B n it is enough to explain 
a method to compute the MGF of the supporting cone at the vertex associated to 
the identity permutation (we denote this by /) and then simply permute the results. 
Nevertheless it is important to stress that, although useful and economical, there is 
no reason to use the same triangulation at each vertex. Similarly, the triangulations 
we use are all regular, but for our purposes there is no need for this property either. 

There are n 2 facets of B n : for any fixed (i,j) : 1 < i, j < n, the collection of 
permutation matrices P satisfying P(i,j) = defines a facet Fi j of B n . Hence, 
every permutation matrix is on exactly n(n — 1) facets and the vertex / is on the 
facets Fi_j,i =/= j. Let C n be the supporting cone at the identity matrix /, then 
the set of facets of C n is {Fij — I}i<i,j< n .i^j- (Note that we need to subtract 
the vertex / from Fjj because the supporting cone is obtained by shifting the 
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supporting polyhedron at the vertex I to the origin.) We are going to apply our 
method CMGF to find the MGF of C„. 

3.1. Step 0: A good projection. C„, as well B n , lie in the Tridimensional space 
K™ = {n x 77. real matrices}. But they lie in different affinc subspaces (the vertex 
of C n is the origin). Let V n be the subspace of W 1 spanned by C n . It is easy to 
see that 

n n 

(3.1) V n = {M G M™ 2 | J2 M ^ k ) = J2 M( - k '^ = 0,V*,i}. 

fc=i fe=i 

Let W n be the vector space R^ 71-1 ) 2 — {(77 — 1) x (n — 1) real matrices}. We define 
a linear map <f> from W 1 to W n by ignoring the entries in the last column and the 
last row of a matrix in M. n , that is, for any M, we define <f>(M) to be the matrix 
(M(i, j))i<ij< n -i- One can check that when restrict </> from V n to W n , 4> is a good 
projection from V n to W n . Let 

C n := 4>(C n ). 

Also, let F,j = 4>(Fi,j) and P = (f>(P), for any permutation matrix P on [n]. 
(These are actually the facets and vertices for A n which is the full-dimensional 
version of B n we explained at the beginning of Section [21) By the linearity of </> the 
facets of C„ are {Pj,j — I}i<i,j<n,i^j, and F^j is defined by the collection of P's 
where P's are permutation matrices (on [n]) satisfying P(i,j) = 0. 

3.2. Step 1: The dual cone D n to C„. The cone C n is full dimensional in 
W = My 1-1 ' . Hence, we can use Lemma [2.11 to find its dual cone. We will first 
define a cone, and then show it is the dual cone to C n . 

Definition 3.1. D n is the cone spanned by rays {Mij}i<i_j< n ^j , where Mij is 
the (n — 1) by (n — 1) matrix such that 

(i) the (i, j)-entry is 1 and all other entries equal zero, if i ^ n and j =/= 77; 

(ii) the entries on the ith row are all —1 and all other entries equal zero, if 
i =/= n and j = n; 

(iii) the entries on the jth column are all —1 and all other entries equal zero, if 
i = 77 and j =/= n. 

Example 3.2 (Example of M^j when 77 = 3). Here we present each 2 by 2 matrix 
Mi_j as a row vector, which is just the first and second row of the matrix in order. 



M lt3 


-1 


-1 








M 2;3 








-1 


-1 


M 3 ,i 


-1 





-1 





M 3 . 2 





-1 





-1 


Ml, 2 





1 








M2.1 








1 






Lemma 3.3. D n is the dual cone to C n inside the vector space W n . 

Proof. For any i, j g [n] and i 7^ j, we need to check that condition (|2.ip is satisfied. 
Note that a ray of C„ is given by the vector P — I, for P a permutation matrix 
adjacent to the identity permutation. Thus it is enough to show that for any 
permutation matrix P on [n], we have (Mij, P) > (Mi t j, I) and the equality holds 
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if and only if P is on F ij; or equivalcntly, P is on the facet F^j. We have the 
following three situations for verification: 

(i) If i ^ n and j ^ n, (Mij, P) is if P is on Fij and is 1 if P is not on Fij. 

(ii) If i y£ n and j = n, (Mi j,P) is —1 if P is on Fi j and is if P is not on 

(iii) If i = n and j ^ n, it is the same as (ii). 

Therefore, D n is the dual cone to C n . □ 

3.3. Step 2: The triangulations of D n . As we mentioned in the last section, 
we will use the idea of toric ideal to find a triangulation of the dual cone D n to 
decompose D n into unimodular cones. 

Lemma 3.4. Let M. be the configuration of vectors {Mij}i<i.j< n ,i^j and [M] 
denote the matrix associated to M, i.e, the rows of[M] are the vectors in A4 written 
as row vectors. The matrix [A4] is totally unimodular, i.e., for any (n— l) 2 linearly 
independent Mij 's, they span a unimodular cone. It follows that all triangulations 
of the cone D n have the same number of maximal dimensional simplices. 

Proof. Up to a rearrangement of rows the matrix [A4] will look as follows: The 
first few rows are the negatives of the vertex-edge incidence matrix of the complete 
bipartite K n —i n —i, then under those rows we have n — 1 cyclically arranged copies 
of an (n — 2) x (n— 2) identity matrix. It is well known that the vertex-edge incidence 
matrix of the complete bipartite K n -\ >n -i is totally unimodular. Moreover it is 
also known, see e.g., Theorem 19.3 in [21], that a matrix A is totally unimodular if 
each collection of columns of A can be split into two parts so that the sum of the 
columns in one part minus the sum of the columns in the other part is a vector with 
entries only 0, +1, and — 1. This characterization of totally unimodular matrices 
is easy to verify in our matrix [M] because whatever partition that works for the 
columns sets of the vertex-edge incidence matrix of the complete bipartite A'„-i,n-i 
works also for the corresponding columns of A4, because the diagonal structure of 
the rows below it. 

The fact that all triangulations have the same number of maximal simplices 
follows from the unimodularity as proved in Corollary 8.9 of [53]. □ 

Therefore, any triangulation of D n gives a decomposition of D n into a set of 
unimodular cones. Since M. defines the vertex figure of D n , it is sufficient to 
triangulate the convex hull of M.. Hence, we consider the toric ideal 

I M ■= (x u+ - x u ~ | uM = 0) 

of M. inside the polynomial ring fc[x] := k[xij : 1 < i, j < n,i j]. Note here 
u G Z™(™ -1 ' is a n(n — 1) dimensional vector indexed by {(i, j) : i,j £ [n], i ^ j}. 

Recall that a circuit of J4 is an irreducible binomial x u+ — x u in 7^ which has 
minimal support. Another result follows immediately from Lcmma l3.3[ Lemma 13.41 
and [221 Proposition 4.11, Proposition 8.11]: 

Lemma 3.5. The set Cm of circuits of the homogeneous toric ideal Im is in fact 
a universal Grobner basis Um f or Im ■ 
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For any partition of [n] = S U T, we denote by us.t G Z™(™ -1 ' the n{n — 1) 
dimensional vector, where 

fl, if i€S,jeT, 
U5,t(»,j) = < -1, if i G T,j G 5, 
[0, otherwise. 

One can easily check that us.t has the following two properties: 

(3.2) u s . T {i,j) + u s ,t(;M) = 0, for any i ^ j. 

(3.3) U5 ^(i,j) + ug r(j, fc) + usj(fc,i) = 0, for any distinct i,j and k. 
We define 

P S ,T ■= - X U ^,T = JJ Xs f - JJ 2^, 

s<ES,t£T seS.teT 

, + .x fl, if*GS,jeT, _ fl, if ieT.jeS, 

where uJ T (i, 7) = < and u s T (i, 7) = < 

I 0, otherwise, I 0, otherwise. 

Proposition 3.6. The set of circuits of Im consists of all the binomials Ps.t's: 

Cm = {Ps.t I S U T is a partition of [n]}. 

Example 3.7. For n = 3, we have 

Cm 



{ P {1}.{2.3} 


= £l,2£l,3 


~ £2,1X3,1, 


P {2,3},{1} 


= X2S%3,1 


- £l,2Xl,3, 


P {2},{1,3} 


= X 2 .lX 2 ,3 


- £1,20:3,2, 


P {1,3},{2} 


= Xi. 2 X 3i2 


- £2,l£2,3, 


P {3},{1,2} 


= X 3 ,lX 3 ,2 


- £l,3£2,3, 


P {1,2},{3} 


= Xi t3 X 2 ,3 


- £3,l£3,2} 



We break the proof of Proposition 13. 61 into several lemmas. Before we state and 
prove the lemmas, we give a formula for the entries in 

i,jS[n\,ijtj 

For any i, j G [n — 1], at most three members of M. are nonzero at (i,j')-entry: 
Mij(i,j) = 1 (this one does not exist if i = j), Mi in (i,n) = —1, and M n j(n,j) = 
— 1. Hence, 

-u(i,n) -u(n,j) i = j; 

u(i, j) - u(i,n) - u(n,j) i ^ j. 
Therefore, we have the following lemma. 

Lemma 3.8. 

u(i, n) + u(n, i) = 0, Vi G [n — 1]; 

u(i, j) — u(i, n) — u(n, j) = 0, Vi ^ j £ [n — 1] . 



(<^)(j,i) = 



U.M = if and only if 



Lemma 3.9. For any partition of [n] = S U T, we /icroe i/iai us,t-M = 0. Hence 
Ps.t is in the toric ideal Im ■ 



Proof. It directly follows from (|3.2p . (|3.3[) . and Lemma 13181 □ 
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Lemma 3.10. For any nonzero u 6 Z n ( n 1 - ) satisfying uA4 = 0, i.e., x u+ — x u G 
Im, t/iere exists a partition of [n] = S U T, so t/iat supp(us : t) C supp(u). 

Proof. Wc first show that there exists t G [n], such that either (t, n) or (n,t) is in 
the support supp(u) of u. Let G supp(u), if either i or j is n, then we are 

done. Otherwise, by Lemma 13.81 we must have either (i, n) or (n, j) in supp(u). 

By Lemma 13.81 again, we conclude that (t, n) G supp(u) if and only if (n, t) G 
supp(u). Let T = {t | (i, n) G supp(u) and/or (n,t) G supp(u)} and S = [n] \ T. 
Both S and T are nonempty. Thus S U T is a partition of [n] . We will show that 
S U T is the partition needed to finish the proof. 

supp(u 5iT ) = {(s,t) \ s e S,t e T} U {(t,s) \ s <E S,t £ T}. Hence, we need 
to show that Vs G S, Vt G T, both (s,i) and (t, s) are in supp(u). If s = n, it 
follows immediately from the definition of T. If s =/= n, (s, n) supp(u) since s G - T. 
Therefore, (uAi)(s,t) = u(s, t) — u(n,t), which implies that (s,t) G supp(u). We 
can similarly show that (t, s) G supp(u) as well. □ 

Lemma 3.11. Let u G satisfying uM — 0, and supp(u) = supp(u5 j t) for 

some partition of [n] = S U T, t/ien 3c G Z sitc/i t/iat u = cu^t- 

Proof. Because ujj = — uy : 5, we can assume that n £ S. Fix to £ ^ an d let 
c := u(n,io), wc will show that u = cus.t- Basically, we need to show that Vs G S 
and G T, u(s, t) = u(n, to) and u(t, s) = — u(n, to). We will show it case by case, 
by using Lemma [3~8l and the facts u(s, n) = u(n, s) = when s ^ n and u(to, t) = 
when t ^ to. 

• If s = n, t = to : u(s, t) = u(n, to) and u(t, s) = u(to, n) = — u(n, to)- 

• If s = n,t ^ to : u(s,t) = u(n,t) = u(to,t) — u(to,n) = u(n,to) and 
u(t, s) = u(t,n) = — u(n,t) = — u(n, to). 

• If s 7^ n, t = to : u(s,t) = u(s,to) = u(s,n) + u(n, to) = u(n, to) and 
u(t, s) = u(t , s) = u(t , n) + u(n, s) = -u(n, t ). 

• If s ^ n,t 7^ to : u(s,t) = u(s,n) + u(n,t) = u(n, to) and u(t, s) = 
u(t, n) + u(n, s) = — u(n, to). 

□ 

Proof of Proposition \3.b\ By Lemma 13. 9( Lemma 13.101 and Lemma 13.111 we know 
that 

Cm C {P s .t I S U T is a partition of [n]}. 

Now we only need to show that for any partition SUT, there does not exist another 
partition S'UT' such that supp(us',T') is strictly contained in suppfug^)- Suppose 
we have such two partitions and let G supp(u5 i x) \supp(us' i T')- Then i and j 
are both in S' or T'. Without loss of generality, we assume they are both in S'. Let 
t G T", then (i, t) and (j, t) are both in the support of us',t>, thus in the support of 
u$,T- But the fact that (i,j) G supp(S', T) indicates that one of i and j is in S and 
the other one is in T. Wherever t is in, we cannot have both (i, t) and (j,t) are in 
the support of us.t- Therefore, we proved that each Ps,t is a circuit. □ 

Corollary 3.12. For any I G [n], 



Gri := {Ps,T \ S UT is a partition of [n] s.t. £ G S} 
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is a Grobner basis of M. with respect to any term order < satisfying xej > Xi.k, for 
any i ^ £. Thus, the set of initial monomials of the elements in Gri are 

Ini(Gri) := { x s _ t \ S LIT is a partition of [n] s.t. £ 6 S}. 

s£S,t£T 

Example 3.13. For n = 3, i = 3 : 

Gr e = { -P{2,3} : {1} = ^2,1^3,1 - 2:1,2^1,3, 

P{1,3},{2} = ^1,2^3,2 ~ ^2,13:2,3, 
P{3},{1,2} = X3,lX3,2 - 351,322,3} 

and 

Ini(Gr e ) = {ar 2 ,ia;3,i) ^1,22=3,2, x 3tl x 3) 2}. 

Recall that Arh(£, n) is the set of all ^-arborescences on [n]. For any T G 
Arb(£, n) , we define the support of T to be supp(T) := { (i, j ) \ i is the parent of j in 
and let M(T) = {Mij | e supp(T)} be the corresponding subset of M defined 
in Lemma I3T41 (Note the support of T actually is the same as the edge set E(T) of 
T. We call it support here to be consistent with the definitions of other supports.) 

Proposition 3.14. For any arborescence T on [n], we define Dt to be the cone 
generated by the rays in the set M \ M(T), i.e., Dt = conc(.M \ M{T)). 

Fix any £ € [n], the term order and Grobner basis described in Corollary \3.1S\ 
give us a triangulation of D n : 

Tri e := (D T \ T G Arb(^,n)}. 

Proof. From the theory of Grobner bases of toric ideals in Theorem l2.9[ the maximal 
simplices are given by the set of transversals, all minimal sets a C {(i,j) | i ^ j € 
[n]} such that a (~l supp(m) ^ 0,Vm G Ini{Grg). Now due to the fact that each 
of the initial monomials are in bijection to the cuts of on the complete graph, the 
transversals are indeed given by all possible arborescences 

{supp(T) I T G Arb(£,7i)}. 

One direction is easy: given any arborescence T on [n] with root £, one sees that 
supp(T) is a transversal. We show the other direction: if given a transversal a, we 
can draw a directed graph G a according to a, i.e., supp(G CT ) = a. We let T be the 
set of all i's such that there does not exist a directed path from I to i. T is empty, 
because otherwise m = Yis^T teT Xs - t ^ Ini(Grg) but a n supp(m) = 0. Therefore, 
for any vertex i, there exists a directed path from £ to i. This implies that there is 
an ^-arborescence as a subgraph of G a . However, by the minimality of a, G a has 
to be this arborescence. 

Finally, from Theorem 12.91 we know that the complement of these transversals 
are precisely the set of simplices of the triangulation. □ 

Example 3.15. For n = 3, i = 3, there are only three trees for K3, thus the three 
3- arborescences Ta,Tb,Tc for Kj, are depicted in Figure^ 



Tr H = { D Ta = conc(7W \ M(T A )) = cone({M 1 , 3 , M 2 , 3 , M 3 ,i, Mi, 2 }) 
D Tb = conc(M\M(T B )) = conc({M h3 , M 2 , 3 , M 3 ,2, Ma.i}) 
D Tc =cone{M\M{T c )) = conc({M 1 , 3 , M 2 , 3 , M li2 , M 2jl })}, 
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Figure 1. 3-arborescences 



where Mij is defined as in Example 

3.4. Step 3: The dual cone to Dt- We have given triangulations Trig of D n . 
By Lemma 13.41 we know this gives a decomposition of D n into a set of unimodular 
cones Dt, one for each arborescence T. Hence we can proceed to find the dual cone 
to each Dt inside W n . 

Recall V n is the subspace spanned by the supporting cone C n at the vertex / 
and can be described by (|3.1|) . We will define cone Ct, for each T € Arb(7, n), in 
the subspace V n , then show Ct ■= </>(Ct) is the dual cone to Dt- 

Definition 3.16. For any directed edge e — (,s, t), (s is pointed to t,) we define the 
weight of e to be the n x n matrix w(e) G M™ , whose (s,t)-entry is 1, (t,t)-entry 
is —1, and all the remaining entries are zero. 

Given T an arborescence on [n] with root £, let v be a vertex of T, then there is 
a unique path from £ to v, we define the weight w T (v) of v with respect to T to be 
the summation of the weights of all the edges on this path. 

Definition 3.17. Let T be an arborescence on [n] with root i. For each directed 
edge e — (i,j) not in T, i.e., e ^ E{T), we define 



W 1 



\s) - w T (t) + w(e). 



More precisely, the entries of W ,e are 

'1, if i ^ j and (i,j) £ cycle(T + e) has the same orientation as e, 

-1, if i ^ j and (i,j) £ cycleiT + e) has the opposite orientation as 

-1, if i = j and i is a vertex in two edges of cycle(T + e) with 

W 1,e (i,j) = ^ both edges having same orientation as e., 

1j if i = j and i is a vertex in two edges of cycle(T + e), with 

both edges having opposite orientation of e., 

t 0, in all other cases. 

where cycle(T + e) denote the unique cycle created by adding e to T 

Let Ct be the cone generated by the set of rays {W T ' e \ e g' E(T)} and Ct be its 
projection under <f> ( the map that ignores the last row and last column of an n x n 
matrix): 

Ct 

Proposition 3.18. 



cone({W/ T ' e 
(1) Each W 



e#E{T)}), C T := 4>{C T ). 
e is in the subspace 



V n = {M e M" 2 I ^M(i,k) 



E 



M(fc,j) = 0,Vi,j}. 
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Hence, Ct is in V n . 
(2) Ct is the dual cone to Dt in the vector space W n = R^™^ 1 ) . 

Proof. (1) We observe that for each row or column of W T ' e , there are cither 
one 1, one —1 and the other entries are zeros or all entries are zeros. 
(2) Ct is the cone generated by the set of rays {4>(W T ' e ) \ e ^ E(T)}, and Dt 
is the cone generated by the set of rays {Mjj | (i, j) ^ supp(T)}. Recall 4> 
is the map that ignores the entries in the last column and the last row of a 

2 

matrix in V n C K" . Hence, we have 

(j)(W T ' e )(k,£) = W T ' e (k,l),Vl <k,£<n-l. 

To check whether Ct is the dual cone to Dt, it is enough to check for 
any directed edge e = (s,t) £ E(T) and any supp(T), we have 

{(f>(W T - e ), Mij) is positive when = (s,t) and is otherwise. In fact, 
we will show that (<f)(W T,e ), M^j) — o~(i,j),( s ,t)- There are three situations. 

• If i ^ n and j + n, then (<j)(W T > e ),Tl' hj '} = c/)(W T > e )(i, j) = W T > e {i,j). 

• Iii = n and j ^ n, then {(f){W T > e ), M i;j ) = ^^(-^(^'^(fc^')) = 
E n k Zl(-W T > e (k,j)) = W T ^(n,j) = W T ^{i,j). _ 

• If % n and j = n, similarly we have (4>(W T ' e ), M^j) = W T ' e (i, j). 
Hence, for every situation {(j)(W T,e ), M iy j) = W T ' e (i, j). However, since the 
only edge in cycle(T + e) not in T is e, W T ' e (i,j) = <%,j),(s,tV 

□ 



Example 3.19. When n = 3,i = 3, as before we will present W T ' e as a row 
vector, which is just the first, second and last row of the matrix in order. For the 
3-arborescence Ta in Figure^ we have four directed edges to be added, the edges 
(1,2), (1,3), (3,1) and (2,3). 

W t a ,(i,2) ;_iio 1 -1 

w t a ,(i,3) : _i o 1 1-10 1-1 

W T A ,(2,3) : o -1 1 1 -1' 

W t a ,(3,i) : o -110 1-10 
Similarly we have edges (1, 3), (2, 1), (2, 3) and (3,2) to be added onto the 3- 



arborescence Tb in 
arborescence Tq- 

W T B ,(1,3) . _ X 

W T B ,(2,1) . _ X 

W T B ,(2,3) . _ x 

W T B ,(3,2) . Y 

^0,(1,2) . _ x 

W T C ,(1,3) . _! 

W Tc -( 2 'V : 

M/T C ,(2,3) . o 



Figure [7] and edges (1, 2), (1, 3), (2, 1) and (2,3) for the 3- 




1 
1 

-1 
1 








-1 
-1 





-1 
-1 



1 



1 

-1 
1 
1 

-1 









1 

-1 



1 
1 



-1 



-1 




-1 



-1 



3.5. Step 4: The multivariate generating function of C n . Because each Dt 
in the triangulation of D n is unimodular, so is the dual cone Ct of Dt- By Lemma 
12.71 we conclude that Ct is unimodular as well and the following proposition: 
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Proposition 3.20. Fixing £ £ [n], the multivariate generating function of C n is 
given by 

( 3 - 4 ) f( c ^)= E II a -n 

TeArb(«,n) e£E(T) y 11 ' 



One observes that Equation (|3.4|) is independent of the choice of I. Thus we have 
the following equality. 

Corollary 3.21. For any £ 1 ,£ 2 £ [n], 

E II (l - rTz vl/T,<! ) = E II fl - y] z »' T ' c )' 

reArb(i!,n)e^E(T) 1 11 ; TeArb(£ 2 ,n) e£E{T) K 11 ' 

4. A RATIONAL FUNCTION FORMULA FOR f(tB n ,z) 

In the last section, we obtained a formula for the multivariate generating function 
of the supporting cone C n of the vertex / of B n . Because of the symmetry of vertices 
of the Birkhoff polytope we can get the MFG of the supporting cone of any other 
vertex of B n . 

Corollary 4.1. The multivariate generating function for the lattice points of the 
supporting cone C n (a) at the vertex a, for a a permutation in S n> is given by 



(4.1) f(C n {a),z) = _ £ n.(I3n^ 



TeArb(£,ri) e<£E(T) 



where W T,e o~ is the matrix obtained from usual matrix multiplication of W T ' e and 
the permutation matrix a. 

Proof of Theorem li.il Note that for any positive integer t, the supporting cone of 
tB n at vertex to~ is still the same supporting cone C n (o~) of B n at the vertex a. 
Then the theorem follows from Corollary 12.41 and Corollarv l4.ll □ 



We conclude this section with an example of Theorem 11.11 for our running ex- 
ample. 

Example 4.2. When n = 3,£ — 3, the three i-arborescences are shown in Figure 
[7J In example \3.1<A we have already calculated W T,e 's. By plugging them in, we 
get the three parts of the products of rational functions contributing to /(C3, z) : 



n l - 



1 1 



tE(T A )^ lZ WTA 'l 1 - 21,222,12a, 1^2,2 1 ~ ^l.S^^l^ij^j^j 



1 - 22,323, 22 2 i2 3j 3 1 - 2 2 ,223,l2 2il 2l3 \ ' 



tE{T B ) 



1 - 21,323,12^^3 j 1 - 21,2^2,1^1,1^2,2 
1 1 



1 - Z2,3Z3,lZl,2ZlAZ 2 lz 3 l 1 _ Zl,l23,2*l 2^3 \ ' 



IS 



J. A. DE LOERA, F. LIU, AND R. YOSHIDA 



n - : 1 1 



1 - Zl, 2Z 3 , ]2i {z 3 2 1 ~ 21,323,1% !%3 



1 1 



1 — z 2,l23,22 2i 2Z 3i i 1 — 2 2, 3^3, 2^3^2,2 

Thus, z 11 f(C 3 , z) equals the sum of the three rational functions multiplied by (z\ x z\ 2 z 
In order to compute the same for other vertices we simply permute the results: 

z ta f{C 3 (a),z) = (^ >CT(1) 4 i(7(2) 4 CT(3) ) x 

( 



1 ~ Zi,<t(2)Z2,<t(1)Zi ,a(l) Z 2 ,1(2) 1 ~ Z l,*(3) Z 3,°(2) Z 2Ml) Z l,l(l) Z 2.l(2) Z 3.l(3) 
1 1 



1 ~ Z2,a(3)Z3,<r(2)Z 2 ,l(2) Z 3 ,1(3) 1 _ Z 2, CT (2) z 3,rr(l) Z 2^(l) Z 3,l(2) 
1 1 

1 - Z l,<y(3) z 3,a(l)Z llj(l) Z 3lj{3) 1 - Zl,a(2)Z2,a(l)Z l ail) Z 2 a ^ 
1 1 



1 ~ Z 2,a(3) z 3,a(l)Zl,a(2)Z 1 ] j{1 )Z 2 l {2) Z 3 ^ {3) 1 - Zl,tr(X)^3,<7(2)^l l0r (2)^3,<T(l) 



1 1 

X 



1 - 2i, ct (2)Z3,<t(1) z i ,a(l) z 3 ,a(2) 1 ~ Z l^(3) z 3, ( r(l) 2: ii(i)^3^(3) 

X , — ). 



1 ~ z 2,<r(l)23, ( r(2)22i(2) 2: 3i(l) 1 ~ Z 2,a(3)Z3,a(2)Z 3 >CT (3)2 2 ,<r(2) 

Finally, the summation of all six z ta f (C 3 (o~) , z) gives J{tB 3l z). 

5. The Coefficients of the Ehrhart polynomial of the Birkhoff 

polytope 

In section 5.2 of (4j, Barvinok and Pommersheim derive a formula for the number 
of lattice points of a given integral convex polytopc P in terms of Todd polynomial 
by residue computation of the MGF of P. When P is an integral polytope, their 
formula explicitly indicates formulas for the coefficients of the Ehrhart polynomial 
e(P, t) of P. Especially, this gives us a formula for the volume vol(P) of P, applying 
it we can get Theorem II. 21 We start this section by briefly recalling related results 

in a. 

Definition 5.1. Consider the function 

d 

" 1 - exp(-r&) 

in d + 1 (complex) variables r and £i,...,£z. The function G is analytic in a 
neighborhood of the origin r = £i = . . . = £<j = and therefore there exists an 
expansion 

+ oo 

G(t;£i, ...,&) = ^V'td^&ll < i < d}), 

where tdj({£j|l < i < ti}) = tdj(£i, • • • , £d) is a homogeneous polynomial of 
degree j, called the j-th Todd polynomial in £i,...,£d. H is well-know that 
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tdj({£i|l < i < d}) is a symmetric polynomial with rational coefficients. See page 
110 in j!8j for more information on Todd polynomials. 

Example 5.2. Here are the first three Todd polynomials when d = 3: 

td 3 (a;i,a;2,a;3) = (1/24) (xi + x 2 + x 3 ) (xix 2 + x 2 x 3 + x 3 xi) , 

td a (a?i,a:a,aJ3) = (l/12)x 2 2 +(l/4 : )x 3 x 1 + (l/12)x 3 2 +(l/12)x 1 2 +{l/4)x2X 3 +(l/A)x 1 X2, 
tdi(xi, X2, x 3 ) = (1/2) Xi + (1/2) x 2 + (1/2) x 3 , and as usual tdo(xi, x 2 , x 3 ) = 1. 

Lemma 5.3. (See Algorithm 5.2 in [4] .) Suppose P C K w is d-dimensional integral 
polytope and the multivariate generating function of P is given by 



where £j = { — 1,1}, cij, &t,i, • • ■ , i>i,d £ i/ie cij's are all vertices (with multiple 
occurrences) of P. and cone(&i ! i, . . . , b^d) is unimodular, for each i. For any choice 
of c G such that (c, bi.j) ^ for each i and j, we have a formula for the number 
of lattice points in P : 

(5.2) \P n Z N \ = ^2 =r d 7 , : £ ^77^ td d _ fe ((c, 6 M ), . . . , <c, M). 



Indeed, if we make the substitution .t^ = exp(rci) Formula (|5.ip can be rewritten 



as 



(5.3) /( P, z) = i.£ e 



' exp((c, Qj) 



(1 - exp((c, hs)) exp((c, ' 

Each fraction is a holomorphic function in a neighborhood of r and the d-th co- 
efficient of its Taylor series is a linear combination of Todd polynomials. Thus its 
d-coefficient of the Taylor series is 

((c,6i,i»'--((c,i»i,d>) kl 

Formula (|5.2p is the result of adding these contributions for each rational fraction 
summand. 

It is clear that if Formula (|5.ip is the MGF of an integral polytope P, then we 
have the MGF of any of its dilations: 

Hence, by using the Lemma l5.3[ we get the Ehrhart polynomial of P. 

Lemma 5.4. Suppose P C M. is d-dimensional integral polytope and the mul- 
tivariate generating function of P (produced by Barvinok's algorithm) is given by 
(|5.ip . For any choice of c € M. N such that (c, bij) ^ for each i and j, the Ehrhart 
polynomial of P is 

(5.6) e(P, *)=£^E ^ 777 - ((c, a*})* td d _ fe «c, 6 M ), . . . , (c, fc M )). 



k=0 



n,-=i( c ^ij> 
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In particular, we get a formula for the volume of P : 
(5-7) v°lQP) = -^ ( j C ' a ' ))d ■ 

Proof. Formula (|5.6p directly follows Lcmma [5 . 3l and our earlier discussion. Formula 
(|5.7[) follows the facts that the leading coefficient of e(P,t) is vol(P) and the 0-th 
Todd polynomial is always the constant 1, which can be shown from the Taylor 
expression of the function G(r, £ 1; . . . , ^) defining the Todd polynomials. □ 



Proof of Corollary \1.SX It follows Lemma 15.41 and Theorem 11.11 □ 
To help our readers we wrote an interactive MAPLE implementation of Formula 



(|1.2p in the case of B 3 . It is available at |http : / / www . math . ucdavis . edu/~delo"e ra/RECENT_WQRK/volBirkhof f ; 

Clearly, it would be desirable to apply a suitable variable substitution of Cjj so 
that the expression of the volume has as few terms as possible (preferably keeping 
size of Cij small), with the hope of speeding up the calculations or even in hope of 
finding a purely combinatorial summation. We leave this challenge to the reader 
and conclude with a variable exchange that gives the volume in just two variables 
(it is possible to leave it as a univariate rational function from the substitution 
Cij = iP). If we set c-j = s l P clearly there will be no cancellations. For example 
for 

1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 
1/24 



the case 


n = 


= 3, the volume of B3 equals. 

( st + s 2 t 2 +s 3 t 3)4 


(st 2 +s 2 t- 


st-s 


2 t 2^ s 2 t 3 +s 3 t 2_ 3 2 t 2_ s 3 t 3)(- 3t 3 + s 3 t 2 +s 2 t _ st _ a 2 t 2_ s 3 t 3)(- s 2 t 2 + s 3 t _ 3 2 t _ s 3 t 2-) + 
( st + s 2 t 3 +s 3 t 2)4 


(st 3 + s 2 t- 


st — s 


2 t3)( s 2p +s 3 t 3_ s 2 t 3_ s 3 t 2^ st 2 +3 3 t 3 +s 2 t _ st _ s 2 t 3_ 3 3 t 2j( s 2 t 3 +s 3 t _ s 2 t _ 3 3 t 3) + 
( st 2 + s 2 t + s 3 t 3)4 


( 3t+3 2 t 2_ 


st 2 - 


«'t) (.2t3+ s 3 t _ s 2 t _ 3 3 t 3 \ L(3 + ,3 t+s 2 t 2_, t 2_ s 2 i _ s 3 i 3) f,3 l+e 3(2 _,2 t 2_ s 3 t ) + 
( s t 2 + s 2 t3 + s 3 t )4 


(st 3 +S 2 t 2 


-st 2 


-s 2 t 3 \[s 2 t + s 3 t 3 -s 2 t 3 -s 3 t)(st + s 3 t 3 + s 2 t 2 -st 2 -s 2 t 3 -s 3 t)(s 2 t 3 + s 3 t 2 -s 2 t 2 -s 3 t 3 ) 
( st 3 + s 2 f + s 3 t 2-)4 


(s t + s 2 t 3_ 


st 3 - 


s 2 t)(s 2 t 2 +s 3 t-s 2 t-s 3 t 2 )(st 2 + s 3 t + s 2 t 3 -st 3 -s 2 t-s 3 t 2 )(s 2 t + s 3 t 3 -s 2 t 3 -s 3 t) + 
{st 3 +3 2 t 2 +s 3 t) i 


(st 2 +s 2 t 3 


-st 3 


- s 2 t 2 )( s 2 t + s 3 t 2_ I ,2 t 2_ 3 3 t )(- st + s 3 t 2 + 3 2 t 3„ st 3^ s 2 t 2_ 3 3 t ^ a 2 t 2 +s 3 t 3_ s 2 t 3_ 3 3 t 2^ 
( st + s 2 t 2 + s 3 t 3)4 


(st 2 +s 2 t- 


st — s 


2 t 2 )(st 3 + s 3 t-st-s 3 t 3 )(s 2 t 3 + s 3 t+st 2 -st-s 2 t 2 -s 3 t 3 )(st + s 3 t 2 -st 2 -s 3 t) + 
( 3t+3 2 t 3 +3 3 t 2-)4 


( 3t 3 + 3 2 t _ 


st — s 


2 t 3)( 3t 2 + 3 3 t _ 3f _ 3 3 t 2^ 3 2 t 2 +3 3 t + 3t 3_ 3t _ 3 2 t 3_ 3 3 t 2^ 3t + 3 3 t 3_ 3t 3_ 3 3 t j + 
( 3t 2 + 3 2 t+3 3 t 3-)4 


(st+s 2 t 2 - 


st 2 - 


3 2 t )( 3t 3 + 3 3 t 2_ 3t 2_ 3 3 t 3)( 3 2 t 3 + 3 3 t 2 + 3t _ 3t 2_ 3 2 t _ 3 3 t 3)( 3t 2 + 3 3 t _ 3t _ 3 3 t 2) + 
(at 2 + s 2 f3 +3 3 t )4 


(st 3 + s 2 f2 


-,t 2 


- s 2 t 3 -)(st + s 3 t 2 -st 2 -s 3 t)(s 2 t + s 3 t 2 + st 3 -,t 2 -s 2 t 3 -s 3 t)(st 2 + s 3 t 3 -st 3 -s 3 t 2 ) + 
( 3t 3 + 3 2 f + 3 3 t 2)4 


Ut+S 2 t 3 - 


st 3 - 


s 2 t )( 3t 2 + 3 3 t 3_ 3t 3_ 3 3 t 2-J( 3 2 t 2 + 3 3 t 3 + 3t _ 3t 3_ 3 2 t _ 3 3 t 2)(- 3t 3 + 3 3 t _ 3t _ 3 3 t 3) + 
( 3t 3 + 3 2 f 2 +3 3 t -)4 


(st 2 +s 2 t 3 


-st 3 


-a 2 t 2 )( st + » 3 t3- s t3- 3 3 t ^ 3 2 t+3 3 t 3 + 3t 2_ 3t 3_ 3 2 t 2_ 3 3 t ^ 3t 3 + 3 3 t 2_ 3t 2_ 3 3 t 3^ + 
( 3t + 3 2 t 2 +3 3 t 3-)4 


(st3 + s 3f_ 


st — s 


3 t 3)( 3 2 t 3 +3 3 t 2_ 3 2 t 2_ 3 3 t 3)( 3f 2 +3 3 t _ 3t _ 3 3 t 2)( 3 2 t + 3 3 t 2_ 3 2 t 2_ 3 3 t ) + 
( 3t + 3 2 t 3 +3 3 t 2^4 


(st 2 + s 3 t- 


st-s 


3 t 2^ 3 2 t 2 +3 3 t 3_ 3 2 t 3_ 3 3 t 2")(- 3t 3 + 3 3 t _ 3t _ 3 3 t 3)(- 3 2 t+3 3 t 3_ 3 2 t 3_ 3 3 t ) + 
( 3t 2 +3 2 t+3 3 t 3-)4 


(st 3 + s 3 t 2 


-St 2 


- 3 3 t 3)( 3 2 t 3 + 3 3 t _ 3 2 t _ 3 3 t 3)( 3t + 3 3 t 2_ 3t 2_ 3 3 t )( 3 2 t 2 + 3 3 t _ 3 2 t _ 3 3 t 2) + 
( s t 2 + s 2 t3 + 3 3 f ) 4 


(st + s 3 t 2 - 


st 2 - 


3 3 t )( 3 2 t + 3 3 t 3_ 3 2 t 3_ 3 3 t ^ 3t 3 + 3 3 t 2_ 3t 2_ 3 3 t 3^ 3 2 t 2 + 3 3 t 3_ 3 2 t 3_ 3 3 t 2j + 



6. Integration of polynomials and volumes of faces of the Birkhoff 

polytope 



In this final section we look at two more applications of Theorem going 
beyond the computation of Ehrhart coefficients. 

The first application is to the integration of polynomials over B n . The main 
observation is that, once we know a unimodular cone decomposition for the sup- 
porting cones at all vertices of B n , a formula for the integral (see Formula (|6.ip ) 
follows from Brion's theorem on polyhedra [3J [5] . 

Theorem 6.1. Suppose P C M. N is a d- dimensional integral polytope and the 
multivariate generating function of P is given precisely by Formula (|5.ip . i.e., we 
have full knowledge of a unimodular cone decomposition for each of P 's supporting 
cones and its rays bij and vertices a% . Then for any choice of y 6 M. N such that 
{y, bij) ^ for each i and j, we get a formula for the integral the p-th power of a 
linear form over P 

(6.D / (y, x y dx = 7 J- 1 ?" , ^E^# ai))P+d 



P 



(p+l)(p + 2)...(p + d)^ ■ J\ d (yM. 3 ) 



Notice that although each term in the sum has poles, the poles cancel and the 
sum is an analytic function of y. Since the pth-powcrs of linear forms generate 
the whole vector spaces of polynomials one obtains, from Thcorcm ll.il a formula 
of integration for the polynomial functions over B n (or for that matter, for any 
integral polytope for which we understand their cone decomposition). 

The next application to the computation of Ehrhart polynomials of faces of 
B n . One can easily obtain from Theorem 1 1 . 1 1 similar formulas for the nonnegative 
integral semi-magic squares with structural zeros or forbidden entries (i.e. fixed 
entries are equal zero). Note that any face F of B n , being the intersection of 
finitely many facets, is determine uniquely by the set of entries forced to take 
the value zero. To obtain a generating function for the dilations of a face F of 
B n , f(tF,z), we start from our formula for f(tB n ,z) in Theorem 11.11 For those 
variables xij mandated to be zero, we select a vector A with entries A,j > so 
that the substitution Xij := s Xij docs not create a singularity (this A exists e.g., 
by taking random values from the positive orthant). Call g(tF,z,s) the result of 
doing this substitution on /(££>„, z). We will eventually set s = 0, but first, let 
us check that will not give any singularities. On the numerator s can only appear 
with a nonnegative exponent. It can potentially create a singularity if it appears 
in a factor of the denominator with negative exponent. But, if this occurs, s can 
be factored out and put with a positive exponent in the numerator. Thus, we can 
safely resolve the singularity. Now, we set s = in g(tF, z, s), those terms that had 
a power of s in the numerator disappear. We obtain a multivariate sum of rational 
functions that gives us only the desired lattice points inside tF. We have now two 
examples of this method. First we apply it to obtain a table with the Ehrhart 
polynomials for a (any) facet of B 3 , B4, B 5 , B 6 . 
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n 


Ehrhart polynomial of a facet of B n 


3 


1 + + t 2 + 


4 


1 _l_ 471 + _l_ 15 ^ 4 f 2 _|_ 73 f 3 , 161 ,4 , 83,5 , 61 ,6 , 1 ,7 _i_ 11 ,8 
140 "i" 315 16 1 ~r 60 80 ~t~ 240 28 5040 


5 


1 i 1752847 , , 904325 ,2 , 147579347 f 3 , 8635681 ,4 , 6412937357 , 5 , 18455639 ,6 , 
^ 360360 " 1 77616 L 1 8072064 1 1 415800 L 1 359251200 1 1555200 ~r 

1611167963 ,7 , 95702009 ,8 , 365214839 f 9 , 5561 ,10 , 52388227 ,11 , 42397 ,12, 
261273600 1 38102400 1 457228800 1 28350 1 1437004800 1 8553600 " + " 

4342517 ,13 , 22531 ,14 . 188723 j.15 
9340531200 1 838252800 1 261534873600 


6 


1 i 87450005 , , 102959133218657 ,2 , 14843359199161 ,3 , 2 30f :.i20! >2 8072 8 8 2 1! H > , 1 , 2:57290 -185580 150129 4.5 , 
± 13728792 1 1 4917307 185120 1 322075353600 1 1 3011404556160000 1 1 2 305 3 2 1 3 ! H .18 38 100 ^ " + " 

15435462135033037,6 , 108878694347719 4.7 , 43930*2 1***M>57369 ,8 1 17GGG8 1 77359 1 807 j.9 , 
1448 15 595 72-1800 L 1 11640079 10 190 ' 1 0-102 3 73 705 72 8000 ' 1054508610855200 L " + " 

434798171323757,10 , 1047553900202141,11 , 25028 193 1507924171,12 , 28330897394929,13 , 
19527987228800 1 105450861035520 1 66283398305 1 8 1000 1 23176013414400 " l ~ 

2229552439625171 4.14 , 6610306048279,15 , 21593 1508972 151 ,16 , 2045239925737 ,17 
6628339836518400 1 1 8 1300088828 116 L 1 1406011 180 1736000 1 814847562547200 ' " l " 

1729908621731 4-18 , 21042914689 ,19 , G 1 3892 1 5 2 1069 ,20 , 139856666897 ,21 , 
5121898964582400 1 572 1 17531335680 " 1 194632 160(15 11312000 " 1 (18 1 2 1 2 502 2 89 1 5 9200 " r 

47580345877 ,22 , 4394050999 ,23 , 9700106723 ,24 
4995558790122700800 " 1 1 560 788893205 750 1 GOO 1 246209083227-175968O000 



With the same method we have also computed, for the first time, the Ehrhart 
polynomials for the Chan-Robbins-Yen Polytope CRYi,CRY4,CRY5,CRYe, and 
CRY 7 (see [TJ]): 



n 


Ehrhart polynomial of the Chan-Robbins-Ycn polytopcs 


3 


1 + l£t + i 2 + 


4 


1 _i_ 157 + _i_ 949,2 , 4,3 , 13 ,4 , 1 ,5 , 1 ,6 
1 60 360 3 36 "^20 " 1 " 360 


5 


i i 2843,, i 1087,2 , 16951 4-3 , 723869 ,4 , 1927 4-5 , 2599 ,6 , 113 ^7 , 257 ,8 , 1 ,9 , 1 ,10 
~*~ 840 ~*~ 224 ^ 4320 ^ 362880 1 2880 1 17280 1 5040 1 120960 1 8640 1 362880 


6 


1 , 1494803 , , 15027247,2 , 361525133 ,3 , 364801681 ,4 , 45175393 ,5 , 4314659,6 , 4392257 ,7, 
± ~*~ 360360 1 1965600 1 43243200 L 1 59875200 1 14370048 L 1 3628800 1 13063680 " + " 

781271 4-8 | 75619 ,9 , 15257 ,10 , 22483 ,11 i 29 ,12 , 23 ,13 , 1 ,14 , 
10886400 1 6531840 1 10886400 1 179625600 1 3628800 L 1 66718080 1 111196800 " T " 

1 ,15 
9340531200 


7 


-, | 571574671 , , 41 125 188103 ,2 , 88102713015001 ,3 , 26256060764993 ,4 , 433329666631051 ,5, 
116396280 1 3700495200 " 1 58063725 12000 1 1852538688000 1 44-10092 8 5 12000 ~ t ~ 

615428916451 ,6 , 97984316095277 ,7 , 7939938012827 ,8 , 66150911095 291 4-9, 
120708403200 " 1 470702772 J.8000 L 1 1 1 7690693 1 2000 1 1 370010217984000 

71471423831 ,10 , 4077796979 ,11 , 8513133061 ,12 , 468626303 +13 , 26270857 ,14 , 
1931334451200 1 643778150400 1 965G67225GO00 1 4707627724800" 1 289 700 1 7081)0 " r 

124270847 ,15 , 2371609 ,16 , 1182547 ,17 , 593 ,18 , 149789 ,19 , 
188305108992000 1 62768369664000 1 71137 1856192000 1 10944228556800 " 1 1210 15100 108832000 

2117 ,20 , 1 ,21 
1216-15100 108882000 1 8688905 7 10 188000 



We conclude with some remarks. First, it is natural to ask whether one can 
derive our volume formula from a perturbation of the Birkhoff polytope and then 
applying Lawrence's formula for simple polytopes |20j . We found such a proof using 
a perturbation suggested by B. Sturmfels, but the proof presented here yields more 
results, for example, Corollaries 13.211 and 14.11 can only be obtained this way. Sec- 
ond, it is well known that Brion's and Lawrence's formulas can be proved from the 
properties of characteristic functions of polyhcdra under polarity (see e.g. Corollary 
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2.8 in g] or Theorem 3.2 [6]). On the other hand P. Filliman [III [19] expressed the 
characteristic function any convex polytope P containing the origin as an alternat- 
ing sum of simplices that share supporting hyperplancs with P. The terms in the 
alternating sum are given by a triangulation of the polar polytope of P. Filliman's 
machinery yields in a limiting case Lawrence's volume formulas. Different choices 
of triangulation of the polar of P yield different volume formulas for P. Using Fil- 
liman's duality G. Kupcrberg found (unpublished) other special volume formulas 
that follow from pulling triangulations of the dual of B n . 
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